function [C]=TsallisCost(mu,p,gamma_g,gamma_b,sigma)

C=(p*(1/(1-sigma))*(1-(gamma_g^sigma)-(1-gamma_g)^sigma)+(1-p)*(1/(1-sigma))*(1-(gamma_b^sigma)-(1-gamma_b)^sigma)-(1/(1-sigma))*(1-(mu^sigma)-(1-mu)^sigma));
end

